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Abstract 

We show that rotating shallow water dynamics possesses an approximate (adiabatic-type) pos- 
itive quadratic invariant, which exists not only at mid-latitudes (where its analogue in the quasi- 
geostrophic equation has been previously investigated), but near the equator as well (where the 
quasigeostrophic equation is inapplicable). Deriving the extra invariant, we find "small denomina- 
tors" of two kinds: (1) due to the triad resonances (as in the case of the quasigeostrophic equation) 
and (2) due to the equatorial limit, when the Rossby radius of deformation becomes infinite. We 
show that the "small denominators" of both kinds can be canceled. The presence of the extra 
invariant can lead to the generation of zonal jets. We find that this tendency should be especially 
pronounced near the equator. Similar invariant occurs in magnetically confined fusion plasmas and 
can lead to the emergence of zonal flows. 

Keywords: Rossby waves; Drift waves; Triad resonance; Shallow water; Conservation; Adiabatic invariants; 
Zonal jets; Zonal flows in plasmas 
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I. INTRODUCTION 



To introduce the topic of present paper, let us start with two different physical situations, 
which are known. 



The 



The first is the inverse cascade of energy in two-dimensional hydrodynamics 
inverse cascade is related to the presence of an additional (compared to the 3D hydrody- 
namics) positive quadratic invariant — enstrophy. 

The second situation |3J is the appearance of longer waves in sea^wave turbulence: The 
length of typical waves on the sea surface is often much bigger than those generated directly 
by the wind and increases with time (a process known as wave aging). The sea waves 
generated directly by wind produce — via nonlinear interaction — longer waves, the latter 
produce even longer waves, and so on. This process is a manifestation of the inverse cascade 
^. As in the first situation, the inverse cascade is related to the presence of another (in 
addition to the energy) positive quadratic invariant, in this case the wave action. The 
later invariant holds because the gravity wave dispersion forbids 3-wave interactions, and 
so the main resonance interaction involves 4 waves and conserves the total wave action. 
The conservation of the wave action is similar to the conservation of the total number 
of molecules in rarefied gas (when the main interactions are binary collisions). However, 
unlike the number of molecules, the wave action is only an approximate invariant, whose 
conservation fails in higher order interactions (e.g., 5- wave interactions are possible, and 
they fail to conserve the wave action). 

The present paper considers wave dynamics in rapidly rotating geophysical fluids. An 



Jets j5|], see also the collection of 
Zonal jets well, observed 



integralpart of this dynamics is the emergence of zona 
papers {g]; the stripes on Jupiter make a famous example tJ. 
in the dynamics of magnetized plasmas (which is mathematically similar to geophysical fluid 
dynamics); they appear to act as transport barriers in tokamaks jsl. Plasma regimes with 
zonal jets have become an integral part of modern controlled nuclear fusion installations, in 
particular, ITER. 

It is interesting to see if the emergence of zonal jets can be related to the existence of an 
additional invariant (similar to the two examples given above). This is the main motivation 
of the present paper. 

It is believed that zonal jets emerge as a result of Rossby wave interactions, see 0, 



lOj. The nonlinear dynamics of Rossby waves in the beta plane is often modeled by the 
quasigeostrophic equation (see e.g. 



{AiP - a^il))t + + i^y^i^x - i^xMy = (I.l) 

for the stream function ip{x, y, t) of the horizontal fluid velocity (u, v) = {ipy, —ipx)- Here, a 
is the inverse Rossby radius of deformation, and f3 is the beta parameter characterizing the 
variation of the Coriolis force. The subscripts x, y, t denote partial derivatives, and A is the 
two-dimensional Laplacian. 

It has been shown that the quasigeostrophic equation (II. ip indeed possesses an approx- 
imate (adiabatic-type) quadratic invariant, which requires the inverse cascade to transfer 



12-ll6| 



energy not just to large scales but specifically towards zonal flow 

However, in several physical situations (including Jupiter) zonal jets are well pronounced 
near the equator, while equation fll.l|) is not applicable there. It is a major challenge to see 
if the approach to zonal jets based on the extra invariant works near the equator as well. 

To account for the equatorial region, we consider the rotating shallow water dynamics in 
the beta-plane (equatorial or mid-latitudinal) 

ut + uu:^ + vuy - f{y)v = -gH^, (I.2a) 
vt + uv^ + vvy + f{y)u = -gHy, (I.2b) 
Ht + (Hu), + iHv)y = 0, (I.2c) 



3.g. HP 



e.g. [iTI, |l8| . This system of equations describes the evolution of the horizontal fluid velocity 
{u,v) and the fluid height H (flat bottom is assumed). The function f{y) is the Coriolis 
parameter, and g is the acceleration due to gravity. 

Considering perturbation expansions (see below) for the system f ll.2|) . we find "small 
denominators" not only related to the resonance triads [like in the case of (II. ip ]. but also 
related to the equatorial limit (/ — i- 0). We show that "small denominators" of both kinds 
can be canceled. The possibility of such cancelation is a remarkable property of the rotating 
shallow water system. We are unaware of any other system with similar attributes, even 
remotely. 

Once we pass from a single equation (II. ip with constant coefficients to a system ( 11. 2p of 
three equations with ?/-dependent coefficients, we also face two other problems: 



1. There could be resonant interactions and energy transfer between the Rossby waves 
and the inertia-gravity waves. 

2. The translational symmetry is broken. 

The first problem is resolved due to a general fact observed in a variety of rotating fluid 
systems: The coupling constant in the triad interaction between a slow mode and two fast 
modes vanishes in the equation for the slow mode [3]; in our case, the slow is the Rossby 
mode, and the fast is the inertia-gravity mode. The second problem makes perturbation 
expansions significantly harder; in particular, the Rossby mode needs to be refined (Section 
IVAD. 

We show that in the limit of weak nonlinearity, the system ( 11.21) possesses an additional 
approximate (adiabatic-type) invariant, which is described in Sec. [TTl Before its formal 
derivation in Sec. llVt we demonstrate how the presence of this invariant makes the inverse 
cascade anisotropic and steers energy toward zonal flow (Sec. IIII Al) . Even more specific 
features, observed in some experiments, are in agreement with the proposed picture: 

• Near the equator, the emergence of zonal jets is more pronounced than it is at mid- 
latitudes fSec HTTBl) . 

• In the opposite limit (when typical length of waves excited by forcing is much greater 
than the Rossby radius of deformation), the extra invariant just says that the energy 
should transfer into the sector of wave vectors k with polar angles > 60° (Sec. IIII CI) . 

During the last half century, several ideas were proposed to explain the emergence of zonal 



jets through the dynamics of weakly nonlinear Rossby waves, e.g., random wave closures 



wave kinetic equation [20|, modulational instability [21], and almost resonant interactions 
22| . Since these approaches consider the weakly nonlinear regime, we believe, they should 
be intimately related to the presence of the extra invariant. For the reasons discussed above, 
it is crucial to see that this invariant is present not only in the quasigeostrophic equation, 
but in the shallow water system as well. 

To explain often powerful equatorial zonal jets, a deep approach was developed over the 
past 40 years; it derives the formation of zonal jets from the instability of the equatorial 



mixed Rossby-gravity waves (see 23|, |2J] and references therein). Since this approach is 



also based on small nonlinearity, the presence of the extra invariant remains relevant. We 



emphasize that, irrespective of the detailed initial instability mechanism, it is the weak 
turbulent inverse cascade that controls the geometry of any emergent large-scale feature. 
Our theory then provides an underlying dynamical mechanism specific to the emergence of 
zonal jets, as observed in numerical simulations. We do mention, however, the following issue 
that merits future investigation. Very near the equator (within a couple of degrees), there is 
evidence that dynamical terms neglected in the "traditional approximation" (in which the 
horizontal component of Earth's angular velocity is ignored, producing the standard beta 
plane approximation; see {25!) may play a significant role 24 1. Our theory ignores such 
terms, so it is presently unknown if the extra invariant still exists under these conditions. 
The present work, at minimum, extends the domain of validity of the extra invariant from 
mid-latitudes to immediate vicinity of the equator, where the quasigeostrophic equation is 
already inapplicable, but the "traditional approximation" still remains valid. 

We should also mention several mechanisms that connect zonal jets with strong nonlin- 
earity, when nonlinear terms (in dynamic equations) are similar in magnitude to linear terms 
(see, e.g., j26l-l28|). Note that situations where linear terms are negligible compared to the 
nonlinear terms are not related to the formation of zonal jets because only the linear terms 
are anisotropic (the phenomenon of spontaneous emergence of anisotropy is not relevant 
here because the emerging jets are always observed to be parallel to the equator). As long 
as the linear terms are significant, even if not dominant, the extra invariant should still play 
an important role. 



II. EXTRA INVARIANT 



It is well known that ( ]I.2|) conserves the following quantities: 



o The energy, defined explicitly below by equation (IIV.SP 



o The infinite set of potential vorticity integrals 

' vx -Uy + f{y) 
H 



H dxdy, 



(II.l) 



where F is an arbitrary function of a single variable. These conservation laws are 
related to the advective conservation of potential vorticity |17l |. 



In particular, for F = 1 one obtains the total mass J H dxdy; its conservation implies 
the existence of time-independent space-averaged fluid height H, such that 



{H - H) dxdy = 0. (II.2) 
o The x-momentum 

J uHdxdy, (II-3) 
which is related to translational symmetry in zonal direction. 
We will see that the dynamics fll.2p adiabatically conserves three more quantities. 



Before we describe them, let us eliminate some dimensional parameters by rewriting (II. 2p 
in terms of the fractional relative height 

h{x,y,t) = (II.4) 

and rescaling 

ct^t, {u/c,v/c) ^ {u,v), f{y)/c-^ f{y), where c= \fgS. (II.5) 

Then the shallow water dynamics (II. 2^ takes the form 

Ut + uu^ + vuy- f{y)v = -h^, (11.6a) 

vt + uv^ + vvy + f{y)u = -hy, (II.6b) 

ht + Ux + Vy + {hu)x + {hv)y = 0, (II. 6c) 

where u, v, h are dimensionless, while x, y, t and 1// have dimension of length. 



Consider the linearized perturbational potential vorticity (see [29|]) 

Q = v,-Uy-f{y)h, (II.7) 
which, according to (111.61) . obeys the equation of motion 

Qt + (uQ). + ivQ)y = -Pv{l + h), where ^ = f'iy). (II.8) 

Since our goal is to describe the energy transfer in Fourier space, we consider the Fourier 
transform Qk of the field Q: 

Q{x, y,t) = J Q^{t) e^(^^+'?^) dpdq [k = (p, q),e=p^ + q% (II.9) 



We will show that the shallow water dynamics (1II.6P adiabatically conserves three quan- 
tities of the form 



I= \ / Xk QkQ-kdpdq. 



(11.10) 



The first is the energy of the Rossby waves (the inertia-gravity component is excluded); it 
has 

1 



-j^energy _ 

The second is the enstrophy of the Rossby component; it has 



X enstrophy -i 

In addition to these two, there is an extra invariant with 
1 



(Il.Ua) 



(Il.llb) 



X extra 
k 



, f{q + pV^) , f{q-pV3) 2V3fp 

arctan — arctan - 



k^ 



P + k^ 



(II.Uc) 



This expression is nonsingular as / — )■ 



extra 



5p^ + 42p2g2 + 21g^ 



f + 0{f). (11.12) 



The notion of adiabatic conservation here is similar to that in the theory of dynamical 
systems 30| : The adiabatic invariants are conserved approximately over long time. However, 



here adiabatic conservation is due not to the slowness of parameter change in time, but to 
the slowness of spatial change and to the smallness of the wave amplitudes. This adiabatic 
conservation is due to the presence of two small parameters (see Section flV BP : First is the 
strength of nonlinearity, compared to the beta effect [see (IIV.IOI) ]. and second is the degree 
of spatial inhomogeneity, i.e., the slowness of the dependence of the Coriolis force on the 
latitudinal coordinate y, compared to the length scale L of field variation and the Rossby 
radius of deformation [see (1IV.11|) ]. 

In the present paper, we derive adiabatic conservation of the above integrals in the low- 
est possible non-trivial — leading — orders. We aim here only to establish the fact of adia- 
batic conservation at minimal accuracy (although the actual conservation accuracy might 
be higher, or the conservation time interval might be longer). Our derivation is formal 
asymptotic, but we take special care that no secular terms appear. 

In the present paper we consider the simplest possible case, making two simplifications: 



1. We use the beta-plane approximation (disregarding complications of spherical geom- 
etry). 

2. We assume the fields {u, v, h) vanish at infinity, i.e., at the periphery of the beta plane. 
The same assumption is made for the quasigeostrophic equation OI.ip when considering 
its invariants. Without this assumption, we need to account for the boundary terms. 
These can be dealt with, but their account leads to heavy mathematical calculations, 
which will not be presented here. 

The central result of our derivation — which allows us to establish the extra invariant 
near the equator — is the possibility to cancel "small denominators" at / — )■ 0; see equations 

firosi) . f iroiji . 

We derive the extra invariant in Sec. llVt but first we demonstrate the connection between 
the invariant and zonal jets. 



III. THE EMERGENCE OF ZONAL JETS 



A. Why the extra invariant impUes the emergence of zonal jets 

The approximate conservation of the energy and enstrophy, contained in the Rossby 
component (see Sec. [Tll) implies the inverse cascade of Rossby wave energy. At the same 
time (as we will see now), the presence of the extra invariant ensures the anisotropy of the 
inverse cascade: The energy is transported not just towards the origin, but specifically to 
the region of the k = (p, g)-plane around the g-axis {\p\ <^ \q\), which corresponds to zonal 
jets. 

Indeed, the extra invariant can be written in the form 

£il dpdq , (III-l) 

where is the Rossby wave energy spectrum, and 0k is the ratio of the extra invariant 
spectral density to the energy spectral density 



arctan — - arctan — - — - 



(III.2) 



see (lll.llc|) and (jll.llap . Figure [T] shows a contour plot of the values of the ratio 0k (on a 



logarithmic scale) vs. k. We pose the following question: Is it possible for the energy from 
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The vafues of hg^Jt^f^ 
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FIG. 1. Contour plot of logxo(0k); the ratio (j)]^ measures how much extra invariant is carried per 
unit energy by an excitation with wave vector k = {p,q); see (IIII.lj )- piI.2|) . The plot spans the 
range 0.1 < p < 10, < q < 10, while / = 1. The values of the ratio 0k at the centers of boxes 
A and B differ by roughly a factor of 7000. Therefore, the only way for the energy to transfer 
towards the origin (via the inverse cascade — be it local or nonlocal) is for it to 'squeeze' around 
the g-axis. 

the region A in Fig. [T]be transferred (via the inverse cascade) into the region 5? The value 
of the ratio in the region B is about 7 x 10^ times greater than its value in the region 
A. So, if the transfer A ^ B did occur, the value / of the extra invariant flllLip would 
significantly increase. The only way for the inverse cascade to transfer the energy towards 
the origin would be to transport the energy (on average) along the level lines of the function 
0k- Thus, the dynamics must display anisotropic "Bose condensation": Spatial anisotropy, 
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The values of log^gi^^i^) 




2 4 6 8 10 

zonal wave number p 



FIG. 2. Same as Fig. [U but for an equatorial region, i.e., for small /; for this particular figure, 
/ = 0.03. Comparison of Fig. [1] and Fig. [2] demonstrates that the energy transfer towards zonal 
jets should be more pronounced near the equator than at mid-latitudes. The values of the ratio 
0k at the centers of boxes A and B now differ by roughly a factor of 30 000. 

which is only weakly broken in small scale dynamics, becomes ever more strongly broken on 
large scales. 

B. Why zonal jets should be more clearly observed near the equator 

The difference between the values of (j) in the regions A and B increases as |/| decreases 
(/ near the equator). Figure [2] shows the values of the ratio 0k when / = 0.03. The 
value of in the center of region B is now about 3 x 10^ times greater than its value in 
region A. Therefore, on the equatorial beta-plane, / ~ 0, the inverse cascade is forced to 
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FIG. 3. The quantity f^cj) is a function of k/f while the polar angle 9 is held fixed. The curves 
correspond to the ten rays with 9 = 0°, 15°, 30°, 45°, 60°, 61°, 76°, 89°, 89.9°, 89.99°. Curves with 
9 < 60° (including 9 = 60°) are almost indistinguishable and marked by one label 9 < 60°. If 
9 > 60°, the curves become horizontal when k <^ f. 



transfer energy even closer to the g-axis. 

The more pronounced formation of zonal jets near the equator can be seen quantitatively 
from Fig. [31 which shows the dependence (on a log-log scale) of the ratio (p vs. the wave 
number k at fixed polar angles 6 (due to symmetries, we need consider polar angles only in 
the range 0° < < 90°). The curves shown in Fig. [3] are steeper for large k/f than for small 
k/f: For large k/f they vary as while for small k/f they vary as k~^ if 6 < 60° and as 
k'^ if 6 > 60°. Therefore, during the inverse cascade, the ratio increases more significantly 
with decreasing if / ^ (near the equator). 

For example, if the energy originated in the region k/f> 20, then the inverse cascade 

11 



must transfer this energy (on average) into the sector 89.9° < 6 < 90°. [Indeed, at 
[k/f > 20 and all 6) is less than /V at {6 = 89.9° and k/f 0).] Such 'tight squeezing' 
of energy around the g-axis hardly can be accounted for by the relative decrease of the 
nonlinearity as / — )■ (which might be expected in some situations). 



C. Long-wavelength limit — polar angle 60° 

Now let us consider the opposite limit where k/ f is small. According to Fig. |3l the inverse 

cascade can now transfer energy anywhere into the sector 

60° < ^ < 90°. (III.3) 



This is exactly the sector that was found [3l|| on the basis of satellite altimeter observations 
of the spectra of very long mid-latitude Rossby waves (with periods of several years). The 
sector (1III.3P is clearly visible in the contour plot of log]^o(0k) for small k/f — see Fig. H) it 
shows the values of the ratio 0k when / = 30. The magnitude of 0k drops sharply when the 
polar angle 6 increases beyond 60°; it is clear that, following any level curve beginning at 



larger k, one may approach the origin only through the sector (]III.3|) . 



We see that if the energy is generated at large scales (much greater than the Rossby 
radius of deformation) then the balance argument, based on the extra invariant, does not 
require the inverse cascade to accumulate energy in zonal flows. This conclusion agrees with 



the investigation [32|, which reported "suppression of the Rhines effect" for large /. 



To conclude this Section, we note that the existence of the extra invariant and the balance 
argument (described in Sees. IIII Al IIII B[ IIII Cp holds for a wide class of wave systems with 
Rossby dispersion law. When the nonlinearity is taken into account, then for some special 



forcing, the energy can still concentrate in zonal flows, even with large /, see [33|. In that 
paper it was also found that, in the short wave case (or near the equator, Sec. IIII Bjl . specially 
arranged forcing can accelerate the formation of zonal jets. 

The presented balance argument for the emergence of zonal jets has the appeal that it is 
based on a (previously unnoticed) conservation law. However, this argument crucially relies 
on the assumption of weak nonlinearity. Whether the nonlinearity is weak depends both 
on the forcing strength and on the location of sources and sinks in Fourier space. Physical 
examples often show that the turbulence is weak in the large-scale part of the inertial range, 
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The values of log^gi^^i^) 




2 4 6 8 10 

zonal wave number p 



FIG. 4. When k <^ f, the extra invariant forces energy to accumulate in the sector 60° < 9 < 90° 
(cf. Fig. [1]); for this particular figure, / = 30. The dashed ray marks polar angle 60°. 



in spite of the fact that the energy spectrum becomes infinite when k — t- 0; e.g., consider 



sea wave turbulence [s]. In the case of geostrophic turbulence, the ratio of the magnitude 
of nonlinear terms to the magnitude of linear terms in the quasigeostrophic equation is the 
Rhines number e = yl/(/3L^). [For simplicity, we consider here the short-wave limit, when 
the Rossby radius of deformation is effectively infinite; more refined estimates will be given 
in Sec. IIVBI ] During the inverse cascade the length scale L increases, while the velocity 
scale A stays roughly constant (determined by the energy), and so, e — > 0. 
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IV. DERIVATION OF THE ADIABATIC INVARIANTS 



In this Section, we demonstrate approximate conservation of the quadratic invariants 
fULTOD with the kernels fULTTD . 



A. Refining the Rossby mode 

Dropping the nonhnear terms in (1II.6P leads to the linearized system 

ut = -fv - hx, (IV.la) 

vt = fu-hy, (IV. lb) 

ht = -u^-Vy, (IV. Ic) 
while the linearized perturbational potential vorticity flll.7|) obeys the equation 

Qt = -/3v. (IV.2) 

Because of the ^/-dependence of coefficients in f llV.ll) . we need to refine the Rossby mode. 
Let us add to Q a correction TZ (to be determined below) that is of higher order with respect 
to the parameter (3, to obtain a new field 

s = Q + n (IV.3) 

such that in the linear approximation (IIV.II) the derivative St will be determined by s alone 
(not by u,v,h taken separately or in any other combination, besides s). Calculations show 
that we need to construct TZ such that 

if-A)n = ^ifu + hy). (IV.4) 

Indeed, if TZ is determined by (]IV.4p then 

if - A)st = if - A)iQt + Ut) = -if - A)i^v) + Pifu + hy)u (IV.5) 

and, according to the dynamics fllV.ll) . 

if - A)st = /3s, + 2/3'vy + I3"v. (IV.6) 

Here the right hand side is [3sx + 0(/3^), so that the non-s terms are indeed pushed to 
higher order. Neglecting the higher order terms, we see that the s-mode has the Rossby 
wave dispersion 
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B. Small parameters 



The energy of the system fill. 61) is 

E = ^Hc^ j[{u^ + v^){l + h) + h^]dxdy. (IV.8) 

In the weakly nonhnear hmit, the integrand of the energy (lIV.Sp reduces to u'^ + v'^ + h'^, 
and we assume u, v, h to have the same magnitude A. 

As mentioned above, we exploit two small parameters: First, the field magnitude A 
should be "small", compared to the beta-effect [see fllV.lOP ]. Second, the Coriolis parameter 
f{y) must be a "slow" function of y, so that f3{y) = f'{y) is "small" [see (lIV.lip ], We first 
define the two non-dimensional small parameters when the field variations are characterized 
by a single length scale L, being the same in x and y directions. 

The magnitude of the linearized potential vorticity is given by 

A 11 
Q oc -, where - = - + /. (IV.9) 

L I Ij 

To measure the degree of nonlinearity we consider the ratio e of the magnitude of the 
nonlinear convective terms in (lII.Sp to the magnitude of the linear term: 

Near the equator (where / ^ -^^), the nonlinearity degree becomes the Rhines number: 
e = A/(5L^. 

Small inhomogeneity means that / changes little over the length scale L. The change 
is A/ ^ (5L. Away from the equator. A/ should be compared to /; near the equator it 
should be compared to L^^ . Thus, to quantify the degree of spatial inhomogeneity we use 
the parameter 

h = (31L, and so, A = eh. (IV.ll) 

In more general situations the length scales in x and y directions can be different (which 
is especially relevant when considering zonal jets). Moreover, the dynamics can be charac- 
terized by a wide range of length scales, and they can change in time (they can easily change 
by an order of magnitude during the inverse cascade). To account for different situations, we 
will just keep track of powers of A and /3 (A — )■ 0, /3 — )■ 0). To maintain the condition that 
the field be small in comparison to the beta-effect, we will assume the existence of a small 
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parameter e, such that A oc e(3. [In general, e and b will have a more complex dependence 
on physical scales than (IIV.10|) and flIV.ll|) .] 



When there is a single length scale L, then the 7^-correction in (1IV.3P is 0{AL^^b) and 
is proportional to f3. However, in a general situation, with many length scales, we can only 
guarantee that TZ oc Indeed, for states almost constant in the zonal direction {d/dx = 0) 
and near the equator (/ = f3y) equation ( ]IV.4p becomes 



(3Vn-^ = mu+^), (IV.12) 
which is reduced by rescaling y = y^y to the form 

exhibiting explicitly the ^/P scale of TZ. 

Since the difference between fields Q and s is small (proportional to y/^), we replace Q 
in f UTTOl) by s: 

I^I' = \j X^s^ s_k dpdq. (IV.14) 

C. Supplementing the quadratic extra invariant with cubic terms 

Our central claim is that the increment A/* = I*{t) — /*(0) remains small over long 
times t. However, this does not necessarily mean that /* is small: can oscillate in time, 
similar to the behavior of adiabatic invariants in the theory of dynamical systems. So, we 



use the approach [3J] and supplement the quadratic integral (]IV.14p with a cubic part 

Jsuppl = J* + /cubic. (IV. 15) 

then require J^^'^ppI to vanish to leading order. The general form of the cubic correction is 

/cubic ^ 1 y [y-«^^„2«3 + y-7^^1^2t;3 + Y.fJ'hhh] rfl23 

+ Yl^,fh,hus + Y,f,-hh2Vs] d,23 + j YZ^'u.v^hd.^s (IV.16) 

with 10 kernels y""", F^^^, . . .. Here and throughout the rest of this paper a subscript 
j stands for the wave vector kj = {pj,qj) {j = 1,2,3); e.g. ui = Uki, likewise, F123 = 

16 



Y(ki, k2, ks) for any kernel Y, and S123 = S(ki + k2 + ka), dus = dki dk2 dk^. In addition, 
a subscript — j will denote — k^, in particular, YL123 = F(— ki, k2, ka). 

The form of the shallow water system allows us to consider the following more simple 
form of the cubic correction 

^cubic _ 1 y ^^^^ [Mi23% + N,23V3 + T,23 h] ^123 + ^ ^ ^-123^1^253^123 (IV.17) 



with only 4 kernels M, N, T, and Y instead of the ten kernels in (]IV.16|) . The general 



form (IIV.16P and the simplified form flIV.17p lead to the same final result. A much longer 



calculation demonstrates that the kernels in ( jIV.161) must be related in such a way that the 
terms may be collected in the form flIV.17p . 

When calculating J^^^ppI^ we will have contributions of different nonlinearity orders 

jsuppi _ J* ^j^g linear terms in the equations (IV. 18a) 

+ /* due to the quadratic terms in the equations (IV. 18b) 

_l_ jcubic linear terms in the equations (IV. 18c) 

_l_ jcubic quadratic terms in the equations (IV.18d) 



We will see that the first contribution f lIV.18a|) vanishes automatically. Our goal is to show 
that it is possible to find the cubic correction — with non-singular (uniformly bounded) 
kernels — such that the next two contributions (]IV.18bp and (]IV.18cP exactly cancel each 
other, implying that, indeed, P'^pp^ is determined by only higher order terms (IIV.18dp . For- 
mally, one can always achieve such cancellation for any wave system, but the corresponding 



kernels will generally be singu 



only for very few systems 



34 



ar. The possibility to escape these singularities takes place 



35|. Significantly, our results demonstrate that the rotating 



shallow water system is among them. 

D. The time derivative P"pp' 

We can always assume the obvious symmetries: 

= -^-k, ^123 = ^213) -^123 = -^213) 2^123 = T213, ^123 = ^213 = ^321- (IV. 19) 

According to the rotating shallow water dynamics (III.6p . (III.8p . along with the definition 
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(lIV.Sp . we have to leading orders 



-iQi si+ {-ipiUs - iqiVs) S2 5_i23 



di 



123 



(IV.20) 



[in accordance with our notations (introduced in Sec. IIV Cp . 5_i23 = 5(— ki + k2 H-ka)]. The 
equation flIV.20p exphcitly displays the contributions summarized in (llV.lSb fec). We will 
determine the kernels X, M, N, T, Y from the requirement that the right hand side of this 
expression vanish, and that they be nonsingular. 

First, the integral J XiQi s^iSidi vanishes automatically since is even, and is 
odd. 



Using (ill. 71) and (IIV. 31) . we substitute 



Us 



ipsvs - fh - S3. 



(IV.21) 



into flIV.20p and collect terms into three groups: those containing (1) ssv, (2) ssh, and (3) 



sss: 



/supply 1 J d,2sS^S2 



V3 



[fq3Mi23 - fp-iNus - klTi23, + i{p3,piXi + P3P2^2 + QsQiXi + 9392^2)5123] 



93 

h 

+ — [w?.<li.Mi2i,pm + i{f + 93)^123 - fPzTi2s + f{PiXi + ^2^2)^123] 

93 

^3 

H [i/A^123 - J53T123 + {PlXi + J92^2)5l23] 

93 



-g / ^123 -Si S2 S3 {Vti +Vt2 + ^^3) C?123 



(IV.22) 
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E. Canceling "small denominators" which are due to the equatorial limit (/ ^ k) 



Equating to zero the coefficients of V3 and produces a system of two linear algebraic 
equations, which we solve for the kernels M and A^: 



123 



/ 

m [{PiP3 + 91^3)^1 + {P2P3 + q2q3)X2] + ipiqiXi + q2X2) 



fiP + H) 



6 



123, 



123 



W3T] 



123 



/ 



+ 



(IV.23a) 



(IV.23b) 



W3 [{piP3 + '?ig3)^i + {P2P3 + ^2^3)^2] + ipjpiXi + P2X2] 



6 



123 



These expressions have an apparent singularity when / — )■ 0, which would invalidate our 
perturbational expansion (e.g., making the cubic correction larger than the main quadratic 
part). However, it is possible to choose the kernel T in such a way as to eliminate the 
singularities in M and A^. Both expressions flIV.23P become non-singular as / — )■ if we 
take 



_ PsipiXi + P2X2) + qsiqiXi + ^2^2) 

-tl23 — £9,12 "123- 

This cancels all terms proportional to 1//, and produces 

.f{qiX, + q2X2) , 



(IV.24a) 



M123 

Ni23 



—I- 



P + H 



123; 



' P + kl 



(IV. 24b) 
(IV.24c) 



The denominators in (1IV.24I) still appear to be singular when / — )■ 0, ^3 — ?■ simultaneously. 
However, due to the presence of the delta functions 5i2z = '^(ki + k2 + ka), the condition 
fcs — J- implies ki + k2 — )■ 0, and so, the expressions piXi + ^2-^2 and qiXi + ^2-^2 in the 
numerators (1IV.24P are linear in p3, when /C3 — )■ 0. Therefore, the numerators in (1IV.24P 
are quadratic in /, ^3,^3, and the expressions flIV.24p are bounded. 
Substituting fHv:24l) into fHv:22|l . we find 



jsuppl _ 

~ 2 



^3^1 - Piq3 y , ^3^2 - ^2^3 Y 
^1 + fo , , 9 — ^2 



P + kl 



P + kl 



Si S2 S3 (5i23 d 



123 



^ j ^123 Si S2 S3 (^^1 + ^^2 + ^3) di2Z 



(IV.25) 
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F. Canceling the triad resonance "small denominators" 



Note the obvious identity 



Pi 


P2 




P2 


P3 




P3 


Pi 




(l2 




(l2 


?3 









when ki + k2 + ks 



(to see, e.g., the first equahty, substitute pi = —p2 — P3, qi 
equation). Because of flIV.26|) . equation flIV.251) reduces to 



-q2 



(IV.26) 
— from the last 



jsuppl _ 

~ 2 



pm - piqz 



(Xi - X2) Si S2 S3 5i23 C?123 



P + kl 

Yl23 Si S2 S; + ^2 + a,) rfi23- (IV.27) 

Symmetrizing the first term on the right hand side (over all permutations of the indices 
1,2,3), and using again the identity f lIV.26p . we see that (]IV.27p vanishes if 
Pi 92 - P2qi 



123 



+ ^2 + ^^3 



6 



123 



X2 — Xi ^ X3 



P + kl 



X2 _|_ -^1 ~ ^3 



P + kl 



(IV.28) 



P + kl 

It is apparent that this expression is singular at the points (ki,k2,k3) satisfying the reso- 
nance relations 



kl + ka + k3 = 0, 

^ki + + = 0, 



(IV.29a) 
(IV.29b) 



unless the expression in square brackets vanishes at these points. We will see now that the 
latter is indeed the case. 

On the resonance manifold ( 1IV.29I) . the bracketed expression in (1IV.28P may be put in 
the form 



^1^2 -^2^1 



{piXi + P2X2 + P3^3 



/3 P1P2P3 

To obtain this, the dispersion relation (1IV.7I) must be used along with the identities 



(IV.30) 



Pi 


P2 




P2 


P3 




Ps 


Pi 


when < 


Pl+P2+P3 = 






















^2 




^2 












^Qi + Q2 + ^3 = 0, 



(IV.31) 

which are similar to (]IV.26|) . Thus, to obtain a non-singular form for Y123, we have to 
require that the function pX^ be conserved in the triad resonance interactions, i.e., that the 
equation 

^1X1+^2X2+^3X3 = (IV.32) 
must hold in all points (ki, k2, k3) of the resonance manifold (IIV.29p . 
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G. Kernel X 



The requirement (1IV.32P is satisfied for tlie following five functions: 

P^k = P, 
PX\, = q, 

pX^ = a = I In [fiq + pV3r + k']-l In [fiq - pV?>f + 

pXk = r^k '= arctan f{q + pVs) /k'^ — arctan f{q — pVs) / 

For the functions ^k and ?7k see their physical meaning remains unclear, let alone their 
possible relation to some continuous symmetries. 

For the function (lIV.33al) . the integral / is the energy of the Rossby component. 

The function (lIV.33bp corresponds to the enstrophy. More precisely, in this case the 
integral / is the zonal (East- West) momentum, which is a linear combination of the energy 
and enstrophy. 



(IV.33a) 
(IV.33b) 
(IV.33c) 
(IV.33d) 
(IV.33e) 



The function ( 1IV.33cl) corresponds to the North-South momentum 



However, this 



choice fails to give a physically meaningful quantity in real (coordinate) space because the 
corresponding function Xk is singular (when p — )■ 0). This singularity means that respective 
invariant in real space 

1=1 J Xir,, r2)s(ri, t)s(r2, t)dr,dr2 (IV.34) 



has kernel X(ri, r2) which does not vanish at large separation ri — see 37| for a detailed 
discussion. 

The function (lIV.33d|) fails to produce an invariant either. This is because ^k is even in 
k, and so, is odd, contradicting the symmetry (jIV.191) . 

Unlike ^k, the function ?7k is odd, and the corresponding kernel Xk determines an extra 
invariant for rotating shallow water dynamics. The previously described kernel ( lll.llcl) is a 
linear combination of the functions ( IIV.33aP and (IIV.33ep . 

The proof of the fact that the functions ^k and ?7k are conserved in triad resonance 
interactions has recently been significantly simplified. The new proof is more straightforward 
and can be accomplished with the aid of symbolic algebra software. Indeed, 

if{q + pV^) + k^ 



+ ^'7k = In -^k where Zk 



if{q-pV3) + k^'' 



(IV.35) 
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(In denotes the principal branch of the complex logarithm, with argument between — tt and 
tt), and the required conservation equation 



implies 



(6 + m) + (6 + iV2) + (^3 + ^^3) = 



Zi Z2 Z3 — 1. 



(1V.36) 



(IV.37) 



Now, using flIV.29ap . substitute = —pi — P2, = —Qi — (I2 into flIV.29bp and flIV.37p . 
These equations may then be reduced to two polynomial equations of degree 5 in pi, gi, p2, '?2- 
It is easy to check (e.g., with Mathematica software) that these two polynomials are 
identical up to a constant factor. It follows immediately that the resonance equations fllV.29p 
imply (1IV.37|) . and hence that 



In Zi + In Z2 + In Z^ = 2iTmi, where m = 0, ±1, ±2, . . . . 



(IV.38) 



Continuity considerations require m = 
Thus, there are three invariants: 



13| . and the conservation (IIV.36I) then follows. 



the energy of the Rossby component [corresponding to flIV.33aP 



the enstrophy [corresponding to (lIV.33bp 



the extra invariant [corresponding to (lIV.33ep 



H. Dropping cubic terms 

The cubic terms I™^^'^ have served their purpose in the proof, and can now be dropped. 



similar to the argument 38|| for the quasigeostrophic equation. To see this, first, note that 
the /3^-terms in (]IV.6P can be neglected over a time interval of length at most of order I3~^ 
with u < 2. We also need to consider time intervals containing many wave periods, and so 
V > 1. For specificity, we choose v = 3/2. 

Considering flIV.20p . we have neglected terms oc A^(3^^'^ oc y4^e/3^/^ (such terms come 
from neglecting 7^-correction in the nonlinear terms of the shallow water equations). There- 
fore, over a time t oc /3~^^'^, the error can accumulate at most up to a total error oc A^e. The 



M, N, and T-corrections in pV.17p have the order oc A^efi] the F-correction in pV.17p 
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has the order / [i oc A^e [the kernel Y is proportional to while the kernels M,N,T 
are 0(1), /3 — t- 0]. So, all cubic corrections are within the total conservation error oc A'^e and 
can be safely dropped. As alluded to earlier, these corrections were needed in the derivation 
only to control oscillatory terms; their amplitude is now seen to be small, but their time 
derivative is large (has lower order). 



V. REMARKS 



A. Unique invariant 

The existence of an extra invariant motivates a natural question: Do there exist other 
invariants in the shallow water system? The answer appears to be "No" , although rigorous 
investigation of this question has not been attempted. To elaborate, if such an invariant did 
exist, then the resonance triad interaction flIV.29|) would seem to have another conserved 
quantity, besides (lIV.33p . The latter, however, is known to be untrue 



. This was estab- 



lished by the connection [40[ between invariants of wave interactions and Web geometry 41|. 
It has not been ruled out, however, that the shallow water system (II. 2|) has several invari- 
ants, which collapse into a single invariant for the quasigeostrophic equation (11. 1^ : though 
this seems unlikely. 

The connection to the Web geometry also shows that the dispersion laws that admit 
extra invariants are extremely rare. We are aware of only one other physical system (besides 
Rossby waves) that possesses extra invariants. This is the generalization of the Korteweg-de 
Vries (KdV) equation for two spatial dimensions 

{ipt + i'lpx + 1pxxx)x = i^yy] (V. 1) 

it has dispersion law 

fi(p,g) = V--. (V.2) 
P 

Equation (IV.ID is integrable via the inverse scattering method and has infinitely many extra 
invariants 35|]. The system ( IV. ip is called the Kadomtsev-Petviashvili equation of the first 



kind (KPl); the Kadomtsev-Petviashvili equation of the second kind (KP2) has a minus 
sign in front of the term on the right of (IV.ip : because of this, triad resonances do not exist 
at all for KP2. 
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Unlike to the KPl case, the Rossby wave triad resonance admits only one extra invariant; 



38|. So, 



and moreover, it is impossible to extend this invariant to the next nonlinearity order 
the extra invariant of the shallow water dynamics (II.2p is an attribute of weak nonlinearity. 
The triad resonances that admit finite number of extra invariants are even more rare 



than the ones with infinitely many invariants (see 39|]): The former constitute a several 
parameter family among all functions depending on two variables; and moreover, most of 
the members of this family are not even elementary functions and hardly can be dispersion 
laws of physical systems. 

B. The impact of the extra invariant on statistical equilibrium 

The existence of the extra invariant may provide barriers to statistical equilibration. The 



equilibrium theory for the quasigeostrophic (|42| and references therein) and shallow water 
systems 43| were derived by enforcing only the exact conservation laws (energy, momentum, 
and the potential vorticity hierarchy). Since the latter fully define the equilibrium state 
(under the ergodic hypothesis), the adiabatic conservation laws will generally be violated. 
Given that true equilibration is an infinite time property, the presence of an adiabatic 
invariant does not lead to any mathematical contradiction here. However, there are practical 
issues since the extra invariant could greatly increase the equilibration time scale. This issue 
needs to be investigated. 

We should also note a parallel between the existence of the extra invariant, determined 
only by the Rossby component, and the equilibrium theory. In the latter it is found that 
although the inertia-gravity waves do remove some of the initial energy to small scale surface 
ripples, they do not inhibit the inverse cascade of the remaining energy to form large-scale 
vortex equilibria. 



C. Using perturbational potential vorticity instead of its linearization 

Since the extra conservation holds only in the weakly nonlinear limit, to the same accuracy 
we are free to write the invariants in terms of the perturbational potential vorticity 

^ _ Vx-Uy + f\y) f{y) 
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instead of the linearized perturbational potential vorticity Q, equation (111.70 . We have 
Q ^ HQ (with the error due to nonhnear terms), and instead of OILlOp . 

/=^| X^Q^Q^^dpdq, (V.4) 

where Qk is the Fourier transform of the field Q. 



D. Can rotating shallow water dynamics be approximated by a single equation? 

There is a question whether the shallow water system fll.2p can be approximated near 
the equator by a single equation. Certainly, in the rigid lid approximation {H = const) the 
system (\l.2\i is reduced to the equation of 2D hydrodynamics with beta-effect. However, 
the shallow water dynamics contain three independent variables u, v, H, and accordingly 
the system (II. 2 p contains 3 time derivatives. We allow significant deviations of H from its 
average value H. 

We have attempted to approximate the equatorial shallow water dynamics by a single 
equation 

a, = n,a, + /vi'_„,a.a3i-.3 43 (V.5) 

for the Fourier transform ak(t) of the stream function or some other variable (the notation 
is defined in Sec. IIVCI) . However, we found that such an equation would have insufficient 
accuracy to establish the extra conservation. More specifically, the formula for the kernel W 
would lack one more cancellation in equations similar to (IIV.24I) [numerator in W, instead 
of being quadratic, would be linear in f,p3, qs], and so, the kernel W would be singular. 



E. Possible fast dependence on the y-coordinate 



The extra conservation holds if the coefficients in the shallow water system ( 111.6^ addition 



ally contain fast, but small amplitude, dependence on the y-coordinate. Such inhomogeneity 
may be considered at lowest order as a resonant triad interaction between two Rossby waves, 
with dispersion law pV.7p . and one inhomogeneity wave, with zero dispersion law: 

Pi=P2+P3, (V.6a) 

gi = g2 + gs, (V.6b) 

n{pi,q,) = n{p2,q2)+0; (V.6c) 
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here ki = and k2 = (^2,92) are the Rossby wave vectors, and ks = (^3,^3) is the 

inhomogeneity wave vector. 

If translation symmetry is still maintained in the x-coordinate, one has ^3 = 0. For this 
case, one can readily see that an arbitrary function (p{p, q) that is even in q satisfies 

^{Pi, Qi) = ^{P2, ^2) + (V.6d) 



at each point of the resonance manifold flV.6k bc). Indeed, flV.6ap . with ^3 = 0, and fIV.Gcp 



imply pi = P2 and |gi| = \q2\. In particular, the function (III.llcP is even in q, and the 

conservation ( IV.6dp holds for ip = rj. Thus, the function f lll.llcp is conserved in triad 
resonant interactions of Rossby waves with the inhomogeneity waves. 

Actually, the function (III.llcP is conserved in resonant interactions of any order n (n > 3), 
which involve 2 Rossby waves and n — 2 inhomogeneity waves: 

Pi=P2+ P + + ... + , (V.7a) 

n-2 

qi = q2 + qz + ■ ■ ■ + Qn, (V.7b) 

n(pi, gi) = VL{p2, q2) + P + + ... + ; (V.7c) 

n-2 



Indeed (IV.7aP and flV.7cP imply pi = p2 and |gi| = |g2|5 and therefore, 

Apu qi) = ^iP2, q2) + p + + ... + q (V.7d) 

n-2 

for any function q) which is even in q. 

VI. CONCLUSION 

The Rossby waves have been known [12I, [l^ to possess a rare property: Their triad 
resonance admits an extra conserved quantity: 

ki + ka + k3 = 0, 



fi(ki) + ^](k2) + f](k3) = 

where 



77(ki)+r7(k2)+r7(k3) = (VI.l) 



k = (p, q) =p^ + q^), Q{k) = , (VI.2) 

rj{k) = arctan ( f^^—^^j - arctan ( \ . (VI.3) 
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Despite of the implication (IVI.ip , the extra invariant / is actually independent of the energy 



and momentum (enstrophy) because the integrals fill. 101) or fllll.ll) contain time-dependent 
functions Qkit) or ej^^it). [Recall that = ri(k)/p, 0k = ?7(k)/f2(k).] 
In the present paper, we have established two key results: 

• The Rossby wave extra invariant can be extended to the shallow water dynamics in 
spite of the presence of inertia-gravity waves and in spite of the explicit inhomogeneity 
(the y-dependence of the Coriolis parameter /). 

• The shallow water dynamics possesses an extra invariant in the equatorial limit (when 
/ — )■ 0, but the derivative /' stays away from zero). This limit also leads to small 
denominators, but different from those related to the triad resonance. We have shown 
that it is possible to cancel these small denominators. 

We have also found that for weakly nonlinear shallow water dynamics, the presence of 
the extra invariant constrains the inverse cascade energy transfer to be from small scale 
eddies to large scale zonal flow. The results are also in agreement with some more specific 
experimental features: more pronounced zonal jets near the equator, when / — )■ 0, and 
suppression of zonal jets and the 60° polar angle in the energy spectrum when / — t- oo (see 
Sec. imp . We have seen that the formation of zonal jets is a basic phenomenon that can be 
related to the set of invariants of the rotating shallow water dynamics. 

For future work, it would be crucial to see whether the theoretical predictions agree with 
experimental observations quantitatively, and whether the effects of the extra invariant can 
be clearly resolved from other mechanisms in the plethora of zonal jets phenomena. In 
particular, we believe it important to develop our results for the dynamics of magnetized 
plasmas; it would be very interesting to examine the effects of the extra invariant on the 
formation of internal transport barriers in fusion plasmas. 
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